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ABSTRACT 

Observations of high redshift galaxies have revealed a multitude of large clumpy 
rapidly star- forming galaxies. Their formation scenario and their link to present day 
spirals is still unknown. In this Letter we perform adaptive mesh refinement simula- 
tions of disk formation in a cosmological context that are unrivalled in terms of mass 
and spatial resolution. We find that the so called 'chain-galaxies' and 'clump-clusters' 
are a natural outcome of early epochs of enhanced gas accretion from cold dense 
streams as well as tidally and ram-pressured stripped material from minor mergers 
and satellites. Through interaction with the hot halo gas, this freshly accreted cold 
gas settles into a large disk-like system, not necessarily aligned to an older stellar com- 
ponent, that undergoes fragmentation and subsequent star formation, forming large 
clumps in the mass range 10 T — 10 9 M . Galaxy formation is a complex process at this 
important epoch when most of the central baryons are being acquired through a range 
of different mechanisms - we highlight that a rapid mass loading epoch is required to 
fuel the fragmentation taking place in the massive arms in the outskirts of extended 
disks, an accretion mode that occurs naturally in the hierarchical assembly process at 
early epochs. 
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1 INTRODUCTION 

The morphology and star formation properties of high red- 
shift galaxies are very different from present day quies- 
cent spirals and ellipticals. Large clumpy irregular disks 
with kpc-sized star forming clumps as massive as M& ~ 
10 7 - 10 9 M© ar e observed in the Hubble Ultra Deep 
Field (UDF) (e.g. lElmegreen et alJ I2007L l2009h , a popula- 
tion th at is very rare tod ay. 'Chain galaxies', first identi- 
fied by ICowie et al.l (|l995f ), are believed to be high-redshift 
disky galaxies seen edge-on, while 'clump cluster' galax- 
ies a r e their face on counte rparts ([Dalcanton &; Shectmanl 
1996: lElmegreen et alj|2004h . In optically selected samples, 
high redshift galaxies show very high star fo rmation rates 
up to 100 - 2OOM yr _1 (iDaddi et all l2004h and in re- 
cent spectroscopic observations the y appear to be extended, 
though perturbed, rota t ing d isks (|Forster Schreiber et al.l 
120061 : iGenzel et alJl200dl2008h . The origin of these galaxies 
and how they connect and possibly evolve into present day 
spirals is still unknown. Gas rich maj or mergers give rise to 
large, bulge-dominated rotating disks ([Robertson &; Bullockl 
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2008) even though massi ve clumps can form at l arge radii, 
from globular clusters (IBournaud et al. 2008bh to tidal 
dwar f galaxies (|Elmegreen et al.lll993l : iBarnes &; Hernquistl 
1992). However, major mergers are not frequent enough 
(jDekel et al.ll2009af ) and are more likely to be the origin of 
the rare, extremely high star forming, su b-millimeter galax- 
ies (jZheng et al.ll2004l : Ijogee et al.ll2QQ8h . 

Observational evid e nce (lElmegreen &; Elmegreenl 1 2006 ; 
IBournaud et al.1 l2008al : IShapiro et al.l l2Q08h suggests that 
clumps form in gas rich spiral disks rather than during on- 
going mergers, alt hough the latter scenario can not be com- 
plete ly ruled out (iTaniguchi &; S hiova 2 00ll; lOyerzier et~al] 
l2008h . Re cent work bv IBournaud et al.1 (|20Q7h Thereafter 
B07) and lElmegreen et al.l (|2008h has demonstrated that 
internal disk fragmentation can effectively reproduce many 
of the observables of chain and clump clusters galaxies and 
that these different clumpy systems can have the same origin 
but observed at different inclinations. However, the models 
of B07 still rely on idealized, pre-existing very massive gas 
disks, in order to reproduce the massive clumps and can not 
explain an ongoing, steady-state fragmentation scenario. 

How galaxies acquire their baryons is an open ques- 
tion. The classic picture of galaxy formation within the 
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Figure 1. An RGB-image of the gas showing the disk and accretion region at z ~ 3. The image is constructed using i?=temperature, 
G=metals and i?=density. We can clearly distinguish the cold pristine gas streams in blue connecting directly onto the edge of the disk, 
the shock heated gas in red surrounding the disk and metal rich gas in green being stripped from smaller galaxies interacting with the 
halo and streams of gas. The disk and the interacting satellites stand out since they are cold, dense and metal rich. 



cold dark matter (CDM) scenario assumes that the accreted 
gas is shock heated to the virial temperature, cools radia- 
tively and rains down to form an i nner star- forming rotating 
disk. Recen t theo r etical studies (iBirnboim &; Dekell [2 003 ; 
Keres et al.l 120051: [Pekel fe Birnboiml |2006l: lOcvirk et al l 
' Keres et all (20091 : brooks et alJ bOOa iDekel et all 
have demonstrated that accretion of fresh gas via 
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cold infall can in fact be the dominant process for gas ac- 
cretion for halo masses M < 10 116 M©. In these halos, the 
cooling time for gas of temperature T ~ 10 4 K is shorter 
than the timescale of gas compression and shocks are un- 
able to develop. In halos above this mass, cold accretion per- 
sists as gas is supplied by co ld streams penetra t ing through 
hot m assive halos at z > 2 ([Ocvirk et al.ll2008l : IDekel et al.l 
2009a) whilst the classical hot mode of gas accretion dom- 
inates at lower z. Because of insufficient spatial resolution, 
these studies could not follow the evolution of the accreting 
gas and how the cold streams connect to the central galax- 
ies. The purpose of this Letter is to look in detail at the gas 
accretion and disk formation process using state-of-the-art 
numerical simulations. 



2 NUMERICAL SIMULATION 

We use the ad aptive mesh refinement (AMR) code RAMSES 
(Tevssier 2002) to simulate the formation of a massive disk 
galaxy in a cosmological context including dark matter, gas 
and stars. The gas dynamics are calculated using a second- 
order unsplit Godunov method, while collisionless particles 
(including stars) are evolved using the Particle-Mesh tech- 
nique. T he modelling includes r ealistic recipes for star for- 
mation (Rasera & Tevssier 2006), supernova feedback and 
enrichment ([Dubois &; Tevssierl 120081 ). Metals are advected 
as a passive scalar and are incorporat ed self-consist e ntly i n 
the cooling and heating routine, as in lAgertz et al. (2009), 
and we adopt an initial metallicity of Z = 1O _3 Z0 in the 
high-resolution region. The refinement strategy is based on 



a quasi- Lagrangian approach, so that the number of parti- 
cles per cell remains roughly co nstant, avoiding discreteness 
effects fe.g. lRomeo et al1l2Q08h . The computational domain 
is a 40 Mpc cube containing nested AMR grids of particles 
and gas cells down to a Lagrangian region containing dark 
matter particles of mass m p — 2.2 x 1O 5 M . The effective 
resolution of our initial grid is therefore 2048 3 . We then re- 
fine this base grid according to our refinement strategy, so 
that the maximum resolution is Ax ~ 40 pc in physical units 
at all times. 

For ou r initial conditions we take the Via-Lactea II 
simulation (|Diemand et al.ll20Q8h which forms a Milky Way 
sized dark matter halo that accretes most of its mass (M v i r = 
2 x 10 12 M© at z = 0) by redshift z = 2. We evolved the 
entire simulation to z — at a coarser resolution, here we 
report on the high redshift evolution to z = 2 at which 
point it hosts a disk that is massi ve enough to b e com - 
par ed to the observati ons in e.g. iBournaud et al.l |2p08a) 
and lGenzel et al.l {2006). We use standard galaxy formation 
ingredien ts, with a star formatio n efficiency of 2% (as de- 
fined in iRasera &; Tevssierl l2QQ6h , a star formation density 
threshold nu = 4 cm -3 a nd a supernovae mass lo ading fac- 
tor f w = 10 (as defined in lDubois &; Tevssierl l2008). In order 
to prevent artificial fragmentation, we use a pressure floor 
P ~ 3GAx 2 p 2 , so that we satisfy the iTruelove et al.l (|l997h 
criterion at all times. 



3 RESULTS 

Fig. [T] shows a large scale view of the galactic disk at z ~ 3. 
At this time the dark matter halo has reached a mass of 
Mvir ~ 3.5 x 10 11 M©, while the total baryonic mass in 
the disk (disregarding the bulge) is M bar ~ 2.4 x 10 10 M© 
out of which 50% is gas, putting it in a regime where both 
cold flows and stable shocks can exist (jKeres et al. 2005; 
IDekel fe Birnboiml I200& lOcvirk et all l2Q08h . This striking 
image ties together many aspects present in modern the- 
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Figure 2. A time sequence spanning 40 Myr of the projected gas 
density at z ~ 3 in a 18 x 18 kpc 2 region. The box shows the 
formation of a ~ 10 8 Mq clump via gravitational instability. 

ories of galaxy formation and highlights new complexities. 
Cold streams of gas originating in narrow dark matter fil- 
aments, effectively penetrate the halo and transport cold 
metal-poor gas right down to the proto-galactic disk to fuel 
the star forming region. A comparable amount of metal en- 
riched material reaches the disk in a process that has previ- 
ously been unresolved - material that is hydrodynamically 
stripped from accreting satellites, themselves small disky 
systems, through the interaction with the hot halo and fre- 
quent crossings of the cold streams. 

Streams of cold gas flow into the halo on radial trajec- 
tories, eventually forming orderly rotational motion in an 
extended disk. This gas is in approximate pressure equi- 
librium with the hot halo that has a rotational velocity 
of v YO t ~ 30kms -1 close to the virial radius, increasing 
smoothly to v vo t ~ 200 km s -1 to match the rotation at the 
edge of the disk (r ~ 10 kpc). The ram pressure is signifi- 
cant close to the disk, forcing the streams to curve around 
it. At early times, when the interaction region close to the 
disk is tenuous, streams can 'swing' past the proto-disk be- 
fore being decelerated completely. At later times the turbu- 
lent accretion region carries significant mass and infalling 
cold gas quickly decelerates by plowing through it. We de- 
tect compression and radiative shocks that quickly dissipate 
since the cooling times are very short, resulting in a denser 
configuration for the cold gas. The global outcome of these 
interaction is a turbulent gas heavy disk prone to fragmen- 
tation. Fig. [2] shows a time sequence of the complicated and 
asymmetric gas flows around the gas disk at z ~ 3. The 
figure reveals that many of the large scale spiral arms at 
large radii are not waves, but material arms that can sur- 
vive for an orbital time and that these arms are gravitation- 
ally unstable and can fragm ent into clumps. Grav itational 
instability has been used by Elmegr een et al.l (1993) (here- 
after E93) to explain the formation of massive clumps, as 
large as dwarf galaxies, in the tidal tails of merging galax- 
ies. The typical mass of objects that form within the arms 
is Mj ~ cre ff /G 2 H, where E is the surface density of gas 
within the arm and the effective mass- weighted ID velocity 
dispersion is defined as cr 2 ff — c 2 + cr 2 D where c s is the lo- 
cal sound speed. Using the small region highlighted by the 




Figure 3. Density projection of the stars (left-hand panels) 
and gas (right-hand panels) at z ~ 2.7 illustrating the frag- 
mentation process and the formation of large clumps of mass 
~ 10 7 - 1O 9 M . 

grey square in Fig. [2] we have measured E = 60 M pc -2 
and (Jeff — 25kms -1 , giving Mj ~ 2 x 10 8 M© which agrees 
well with the mass of the forming clump. The internal dis- 
persion velocity is roughly equal to the divergent motions 
across the curved gas filament. The typical velocity disper- 
sion across a A ~ 1 kpc patch of the filament will be of 
the order a c± Xv or h/Tlc — 20 kins" 1 , where the orbital ve- 
locity v or h rsj 200 kms -1 and the curvature radius 1Z C of 
the filament equals the radius of the extended disk. This 
value agrees well with the dominating turbulent component 
of cr e ff. As the interaction region grows in mass and devel- 
ops a more symmetric disk- like morphology we also observe 
massive clump formation in fragmenting spiral waves at in- 
termediate radii. The resulting galaxy is shown in Fig. [3] 
at z ~ 2.7, after many large clumps have formed through 
the above mechanisms. We detect 14 clumps with masses 
between M c \ ~ 5 x 10 7 and 1O 9 M0, of which only the two 
smallest did not form in situ but were infalling satellites. 
The three most massive clumps have M& ^ 7-8 x 10 8 M©. 
In total rsj 15 % of the baryons are in clumps. In the inter- 
action region between the disk and the cold streams, the 
typical arm surface density and velocity dispersion can both 
be estimated using mass average quantities within cylindri- 
cal shells. At (r ~ 8 kpc) we measure (E) ~ 20 M pc -2 
and (a) c± 30 km s _1 , giving rise to clump masses as large 
as Mj ~ 10 9 Mq. Even though we satisfy the Truelove cri- 
terium, convergence in the details of the clump properties 
can be influenced by numerical fragmentation and may re- 
quire more cells per Jeans length. In addition, numerical 
diffusion from bulk flows can lead to an underestimation of 
the turbulent velocity dispersion. Quantifying this is beyond 
the scope of this paper. The detected clumps are located in 
the interaction region between the inner disk and the cold 
streams. In our case this region is not aligned with the ini- 
tial galactic disk, giving rise to a misalignment of the clumps 
with respect to the inner galactic disk (see edge-on images 
in Fig. [3])- Although we believe that this misalignment is 
not typical, it is an elegant explanation f or the formation 
of ' ben t' chains, such as the one reported i n iBournaud et all 
(|2008eJ ). Indeed. [Elmegreen fe Elmegreen! (|2006h report that 
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the typical chain galaxy has clumps mostly aligned in the 
midplane, while in some cases, clumps are seen above and 
below the midplane (outer and inner disk misaligned). In our 
case, the misalignment is due to a third cold stream that is 
perpendicular to the main filament seen in Fig. [T] In a sim- 
ilar scenario, this process has also been invoked to ex plain 
the formation of large polar rings ([Maccio et al.|[2006l ). 

The simulated galaxy is sharing man y properties with 
observed chain and clump cluster galaxies ([Elmegreen et all 
2007). Viewed edge-on, the misaligned disk morphology is 
clearly seen and the overall structure resembles a large 
chain-galaxy. Viewed face on the spiral-like structure has 
a similar morphology as clump clusters or clumpy spirals. 
lEfmegreen &; Elmegreenl ([2005) report that UDF clumpy 
galaxies at z ~ 1.5 — 3 have a stellar mass ~ 6 x 10 10 M© 
and a radius ~ 10 kpc, in striking agreement with our sim- 
ulated galaxy. Not only does the cosmological simulation 
reproduce the observed clumpy morphology and global ro- 
tation of these systems but we also find a realistic metallicity 
gradient and star formation rate of 20 M©yr -1 . The inner 
disk has on average solar metallicity, while that in the clump 
forming region is only ~ 1/10 Z©, due to the accretion of 
pristine gas in the cold streams mixing with stripped satel- 
lite gas. This has the important observational consequence 
that these massive clumps might be devoid of dust, making 
them easier to detect. 

To illustrate 'how disks acquire their baryons', we have 
plotted the mass accretion rate in different gas phases mea- 
sured around our simulated galaxy at z = 5,3 and 2 in 
Fig.g] We define the phases as cold diffuse (T < 2 x 10 5 K, 
n < 0.05 cm -3 ), dense (n > 0.05 cm -3 ), hot diffuse (T > 
2 x 10 5 K, n < 0.05cm -3 ) and stripped (Z > 0.01 Z , 
n < 0.05cm -3 ). Indeed, at z = 3 and 5 the mass accretion 
rates in cold streams is very high (M c± 20 Mqji -1 ). A sig- 
nificant amount of baryons are also accreted from stripped 
satellites, although quantifying this amount is difficult in 
Eulerian schemes since this metal rich material can mix with 
the other gas phases that have never been part of the satel- 
lites. After z rsj 2, the hot mode of accretion dominates, 
making large clump format ion at large radii on l y pos sible 
through galaxy mergers, c. arnes Hernauistl(1992) and 
E93. At z ~ 2, the galaxy has a thin and extended spi- 
ral disk component. Although the gas velocity dispersion 
is still rather high in the disk, the Jeans mass in the spi- 
ral arms is on the order of ~ 10 7 M©, closer to the largest 
giant molecular clouds in present day spiral galaxies. The 
corre sponding gas Q g -parameter ([Goldreich Sz Lvnden-Belll 
Il965h is Q g ~ 1.5 — 2 in the star forming region, indicating 
that the disk is marginally stable and the galaxy has reached 
a quiescent phase with no further large clump formation. 

Fig. [5] shows the dark matter mass accretion rate in the 
simulated galaxy, as a function of time. At z = 2, the ac- 
cretion rate is significantly lower than the average, explain- 
ing why the disk has reached this quiescent phase. A global 
analytical approach for understanding hig h-z disk fragmen- 
tation can be applied (De kel et al.ll2009bf h based on simple 
stability arguments and the disk fraction 5 = Md/Mtot(Rd)- 
Here Md is the baryonic mass in the disk and Mtot(-Rd) is 
the total mass within the disk radius Rd. A (5 ^0.25 — 0.5 
should give rise to large clumps involving a few percent of the 
disk mass and 6 ~ 0.3 — 0.35 is predicted for a steady- state 
fragmentation from moderately clumpy streams. The disk in 




Figure 5. Top panel: dark matter accretion history [and its log- 
arithmic derivati ve (bottom panel)] for t he VL-2 halo using the 
HOP halo finder (Eisenstein & Hut 1998). After a period of mod- 
erate mass increase, during several epochs (e.g. z ~ 3.25 and 
~ 1.75) the halo mass and hence the gaseous mass dramati- 
cally increases. The red rings mark specific times discussed in 
the text. The dotted blue line shows the expected averaged gas 
accretion calculated from extended Press-Schechter (EPS) theory 
(|Neistein et al.ll2006h . 



our simulation at z = 5, 3 and 2 has S = 0.47, 0.33 and 0.33 
respectively which is in excellen t agreemen t with t he above 
prediction (see also Fig. 2 in Deke l et "all (|2009bh V Using 
only the gas in Md renders a lower bound of 5 — 0.17, 0.17 
and 0.1. We point out that the stellar fraction increases sig- 
nificantly towards lower redshifts and this 'hotter' compo- 
nent stabilizes the disk at z ~ 2. 



4 CONCLUSIONS 

We have followed the formation and accretion history of a 
Milky Way sized galaxy using state-of-the-art AMR tech- 
niques. Most of the baryons are in an orderly rotating disk by 
a redshift z = 2, but how they attain this equilibrium is very 
complex and the focus of this work. One of the most impor- 
tant points of this paper is that we can answer the question 
in detail of 'how galaxies get their baryons'. Extending re- 
cent work on the impact of cold streams on ga l axy formation 
(iKeres et al.1 120051: [Dekel fe Birnboiml 120061 : lOcvirk et all 
2008; lDekel et al.ll2009ah , we analyze for the first time how 
single phase narrow cold streams and ram pressure stripped 
debris assembles an extended turbulent rotating disk. Com- 
plex gas interactions takes place in an extended accretion 
region in which infalling gas is decelerated through com- 
pression/radiative shocks and from the pressure gradients 
arising from a hot halo component. 

Prior to z ~ 2, the accretion rate of cold gaseous ma- 
terial onto the disk is the highest and we resolve the grav- 
itational instabilities responsible for the formation of many 
very massive clumps within an extended ~ 10 kpc disk. This 
is about two times larger th an the theoretic al expectations 
of disk sizes at this epoch (|Mo et al.l Il998h . The observed 
morphology, star forming rate, global rotation and metal- 
licity of the system is in good agreement with the observed 
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Figure 4. Mass accretion averaged within spherical shells at redshifts z ~ 5, 3 and 2. The radii are in physical kpc. The lines show the 
total mass flow (solid black) in each shell, cold diffuse (blue solid), hot diffuse (red solid), dense (dotted) and stripped gas (green) (see 
text for definition). We observe a decrease in the overall inflow of material and a change from cold to hot accretion over time. 



clump-cluster and chain gal a xies ([E lmegreen & Elme green! 
l200d : lElmegreen et al.1 120071 : iBournaud et alJ l2008ah . This 
scenario is an extensio n of the disk fragmen tation scenario 
proposed by B07 and lElmegreen et al.l (|2008l ). although here 
studied more consistently within the current cosmological 
framework, and more specifically related to cosmological ac- 
cretion. Therefore, clumpy galaxies should be most frequent 
at this epoch since massive clump formation stops during the 
remaining slow accretion phase and the disk evolves quies- 
cently until z — which will be reported on in a forthcoming 
paper. 
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